Comparison between cell-centered and nodal-based discretization schemes for linear elasticity     

Halvor M. Nilsen  Jan Nordbotten  Xavier Raynaud 《Computational Geosciences》2018,22(1):233-260

In this paper, we study newly developed methods for linear elasticity on polyhedral meshes. Our emphasis is on applications of the methods to geological models. Models of subsurface, and in particular sedimentary rocks, naturally lead to general polyhedral meshes. Numerical methods which can directly handle such representation are highly desirable. Many of the numerical challenges in simulation of subsurface applications come from the lack of robustness and accuracy of numerical methods in the case of highly distorted grids. In this paper, we investigate and compare the Multi-Point Stress Approximation (MPSA) and the Virtual Element Method (VEM) with regard to grid features that are frequently seen in geological models and likely to lead to a lack of accuracy of the methods. In particular, we look at how the methods perform near the incompressible limit. This work shows that both methods are promising for flexible modeling of subsurface mechanics.  相似文献   

Invasions in Marine Communities: Contrasting Species Richness and Community Composition Across Habitats and Salinity     

H. Jimenez  E. Keppel  A. L. Chang  G. M. Ruiz 《Estuaries and Coasts》2018,41(2):484-494

While many studies of non-native species have examined either soft-bottom or hard-bottom marine communities, including artificial structures at docks and marinas, formal comparisons across these habitat types are rare. The number of non-indigenous species (NIS) may differ among habitats, due to differences in species delivery (trade history) and susceptibility to invasions. In this study, we quantitatively compared NIS to native species richness and distribution and examined community similarity across hard-bottom and soft-sediment habitats in San Francisco Bay, California (USA). Benthic invertebrates were sampled using settlement panels (hard-bottom habitats) and sediment grabs (soft-bottom habitats) in 13 paired sites, including eight in higher salinity areas and five in lower salinity areas during 2 years. Mean NIS richness was greatest in hard-bottom habitat at high salinity, being significantly higher than each (a) native species at high salinity and (b) NIS richness at low salinity. In contrast, mean NIS richness in soft-bottom communities was not significantly different from native species richness in either high- or low-salinity waters, nor was there a difference in NIS richness between salinities. For hard-bottom communities, NIS represented an average of 79% of total species richness per sample at high salinity and 78% at low salinity, whereas the comparable values for soft bottom were 46 and 60%, respectively. On average, NIS occurred at a significantly higher frequency (percent of samples) than native species for hard-bottom habitats at both salinities, but this was not the case for soft-bottom habitats. Finally, NIS contributed significantly to the existing community structure (dissimilarity) across habitat types and salinities. Our results show that NIS richness and occurrence frequency is highest in hard-bottom and high-salinity habitat for this Bay but also that NIS contribute strongly to species richness and community structure across each habitat evaluated.  相似文献   

Timing of Proterozoic magmatism in the Sunsas belt,Bolivian Precambrian Shield,SW Amazonian Craton     

Ingrid M.Nedel  Reinhardt A.Fuck  Amarildo S.Ruiz  Gerardo R.Matos-Salinas  Alanielson da C.D.Ferreira 《地学前缘(英文版)》2021,12(6):277-293

We present new U-Pb zircon and monazite ages from the Sunsas belt granitic magmatism in Bolivia,SW Amazonian Craton.The geochronological results revealed four major magmatic events recorded along the Sunsas belt domains.The older igneous event formed a granitic basement coeval to the Rio Apa Terrane(1.95-1.85 Ga)in the southern domain.The second magmatic episode is represented by 1.68 Ga granites associated to the Paraguá Terrane(1.69-1.66 Ga)in the northern domain.The 1.37-1.34 Ga granites related to San Ignacio orogeny represent the third and more pervasive magmatic event,recorded throughout the Sunsas belt.Moreover,magmatic ages of～1.42 Ga revealed that the granitogenesis asso-ciated to the Santa Helena orogeny also affected the Sunsas belt,indicating that it was not restricted to the Jauru Terrane.Lastly,the 1.10-1.04 Ga youngest magmatism was developed during the Sunsas oro-geny and represents the final magmatic evolution related to Rodinia assembly.Likewise,the 1.95-1.85 and 1.68 Ga inherited zircon cores obtained in the～1.3 Ga and 1.0 Ga granite samples suggest strong par-tial melting of the Paleoproterozoic sources.The 1079±14 Ma and 1018±6 Ma monazite crystallization ages can be correlated to the collisional tectono-thermal event of the Sunsas orogeny,associated to reac-tions of medium-to high-grade metamorphism.Thus,the Sunsas belt was built by heterogeneous 1.95-1.85 Ga and 1.68 Ga crustal fragments that were reworked at 1.37-1.34 Ga and 1.10-1.04 Ga related to orogenic collages.Furthermore,the 1.01 Ga monazite age suggests that granites previously dated by zir-con can bear evidence of a younger thermal history.Therefore,the geochronological evolution of the Sunsas belt may have been more complex than previously thought.  相似文献   

Slow quasigeostrophic unstable modes of a lens vortex in a continuously stratified flow     

Hai Yen Nguyen  Richard Schopp  Xavier Carton 《地球物理与天体物理流体动力学》2013,107(3):305-319

This work addresses the linear dynamics underlying the formation of density interfaces at the periphery of energetic vortices, well outside the vortex core, both in the radial and axial directions. We compute numerically the unstable modes of an anticyclonic Gaussian vortex lens in a continuously stratified rotating fluid. The most unstable mode is a slow mode, associated with a critical layer instability located at the vortex periphery. Although the most unstable disturbance has a characteristic vertical scale which is comparable to the vortex height, interestingly, the critical levels of the successively fastest growing modes are closely spaced at intervals along the axial direction that are much smaller than the vortex height.  相似文献   

Wave climate variability in the North-East Atlantic Ocean over the last six decades     

Guillaume Dodet  Xavier Bertin  Rui Taborda 《Ocean Modelling》2010,31(3-4):120-131

Ocean surface gravity waves play a major role in many engineering and environmental problems, both in the open ocean and in coastal zones. Therefore, it is essential to improve our knowledge on spatial and temporal variability of wave climate. This study aims at investigating this variability in the North-East Atlantic Ocean (25°W–0°W and 30°N–60° N), using a 57-year hindcast (1953–2009) obtained with a spectral wave model forced with reanalysis wind fields. The hindcast analysis reveals firstly strong seasonal fluctuations of wave climate, with winters characterized by large and long-period waves of mean direction spreading from south-west to north-west, and summers characterized by smaller and shorter-period waves originating from norther directions. From northern (55°N) to southern (35°N) latitudes, the significant wave height (Hs) decreases by roughly 40%, the mean wave direction (Mwd) rotates clockwise by about 25% while the peak period (Tp) only grows by 5%. These three parameters also exhibit a strong inter-annual variability, particularly when winter-means (from 1st of December to 1st of April) are considered. Linear trend analysis over the studied period shows spatially variable long-term trends, with a significant increase of Hs (up to 0.02 m yr^?1) and a counterclockwise shift of Mwd (up to ?0.1° yr^?1) at northern latitude, contrasting with a fairly constant trend for Hs and a clockwise shift of Mwd (up to +0.15° yr^?1) at southern latitudes. Long-term trends of Tp are less significant, with still a slight increase in the north-eastern part of the study area (up to +0.01 s yr^?1). Eventually, a comparison between the inter-annual variability of the winter-means of the three selected wave parameters and the North Atlantic Oscillation (NAO) reveals: (1) a strong positive correlation between Hs and the NAO index at northern latitudes (correlation coefficient up to R = 0.91) and a significant negative correlation at southern latitudes (up to R = ?0.6); (2) no significant correlation for Mwd north of 40°N and a clear positive correlation southward of 40°N (up to R = 0.8) and (3) a northward increasing correlation for Tp (up to R = 0.8). Long-term trends for Hs, Mwd and Tp are finally explained by a significant increase in the NAO index over the studied period.  相似文献   

Roof-rock contamination of magma along the top of the reservoir for the Bishop Tuff     

Wendell A. Duffield  Joaquin Ruiz  James D. Webster 《Journal of Volcanology and Geothermal Research》1995,69(3-4)

The Bishop Tuff, a well known Quaternary high-silica rhyolite in east-central California, is widely considered the type example of a vertically and monotonically zoned pyroclastic deposit that represents zoning in the source magma reservoir, inverted during the process of pyroclastic emplacement. However, the deposit of plinian pumice, which forms the base of the Bishop Tuff and represents the initial 10% or so of all magma erupted during the event that produced the Bishop Tuff, contains features at odds with monotonie zoning for the reservoir. Relative to overlying ignimbrite, the plinian deposit contains a reversal in trace-element zoning. Moreover, the ⁸⁷Sr/⁸⁶Sr is significantly higher than that in overlying ignimbrite (about 0.7084 vs 0.7064), and melt inclusions trapped in quartz phenocrysts exhibit notable variability of trace-element concentrations, even within a single host crystal (e.g., U: 10.77 to 8.91 ppm).These data have been previously interpreted as due to processes of chemical fractionation and evolution operating within a magma system closed to chemical interactions with its roof rocks. For example, the reversal in trace-element zoning has been explained by the first-erupted magma being erupted from somewhat below the top of a monotonically zoned reservoir. However, we submit that the reversed zoning and other above-noted features can be explained equally well as consequences of minor assimilation of roof rocks into a magma reservoir that was erupted from the top down.The basal part of the Bishop Tuff exhibits extreme concentrations and depletions of trace elements, relative to the average composition of crustal rocks. For example, the upward decrease of Sr in the Bishop magma reservoir (downward decrease in the ignimbrite) results in concentrations as low as 2–4 ppm. Because of the attendant ‘chemical leverage’, assimilation of < 1 wt.% of Sierra Nevada batholith rocks typical of the area could readily reverse an ‘uncontaminated’ Sr (and other trace elements) trend of zoning and could also substantially raise ⁸⁷Sr/⁸⁶Sr. Small-scale trace-element variability in the uppermost part of the Bishop magma reservoir, as recorded by the above-mentioned melt inclusions, may simply reflect melt heterogeneity produced by the process of assimilation.  相似文献   

Hydrodynamical Modeling Of Oceanic Vortices   总被引：1，自引：0，他引：1  

Xavier Carton 《Surveys in Geophysics》2001,22(3):179-263

Mesoscale coherent vortices are numerous in the ocean.Though they possess various structures in temperature and salinity,they are all long-lived, fairly intense and mostly circular. Thephysical variable which best describes the rotation and the density anomaly associated with coherent vortices is potential vorticity. It is diagnostically related to velocity and pressure, when the vortex is stationary. Stationary vortices can be monopolar (circular or elliptical) or multipolar; their stability analysis shows thattransitions between the various stationary shapes are possible when they become unstable. But stable vortices can also undergo unsteady evolutions when perturbed by environmental effects, likelarge-scale shear or strain fields, -effect or topography. Changes in vortex shapes can also result from vortex interactions. such as the pairing, merger or vertical alignment of two vortices, which depend on their relative polarities and depths. Such interactions transfer energy and enstrophy between scales, and are essential in two-dimensional and in geostrophic turbulence. Finally, in relation with the observations, we describe a few mechanisms of vortex generation.  相似文献   

Trace metals in estuarine sediments from the southwestern Spanish coast   总被引：7，自引：0，他引：7  

Ruiz F 《Marine pollution bulletin》2001,42(6):482-490

The impact of river-transported metal pollution and industrial wastes on the metal distribution (Cr, Cu, Zn, Pb) in estuarine sediments was studied in the southwestern Spanish estuaries. Intertidal and subtidal surface sediments of the Tinto-Odiel Estuary are very highly polluted by heavy metals, with geoaccumulation indices up to 4 in the three sedimentary environments studied (channel, channel border and salt marsh). The single exception is the Punta Umbria channel, very protected from the point sources by salt marsh deposits and hydraulic processes. In the remaining two estuaries, pollution (Pb, Cu) was only significant near the harbour situated in the Piedras river mouth, whereas very low values were found in the Guadiana Estuary. In these last rivers, the enrichment factor increases from the channel to the salt marsh sediments.  相似文献   

Evaluation of a high-resolution regional climate simulation over Greenland     

Filip Lefebre  Xavier Fettweis  Hubert Gallée  Jean-Pascal Van Ypersele  Philippe Marbaix  Wouter Greuell  Pierluigi Calanca 《Climate Dynamics》2005,25(1):99-116

A simulation of the 1991 summer has been performed over south Greenland with a coupled atmosphere–snow regional climate model (RCM) forced by the ECMWF re-analysis. The simulation is evaluated with in-situ coastal and ice-sheet atmospheric and glaciological observations. Modelled air temperature, specific humidity, wind speed and radiative fluxes are in good agreement with the available observations, although uncertainties in the radiative transfer scheme need further investigation to improve the model’s performance. In the sub-surface snow-ice model, surface albedo is calculated from the simulated snow grain shape and size, snow depth, meltwater accumulation, cloudiness and ice albedo. The use of snow metamorphism processes allows a realistic modelling of the temporal variations in the surface albedo during both melting periods and accumulation events. Concerning the surface albedo, the main finding is that an accurate albedo simulation during the melting season strongly depends on a proper initialization of the surface conditions which mainly result from winter accumulation processes. Furthermore, in a sensitivity experiment with a constant 0.8 albedo over the whole ice sheet, the average amount of melt decreased by more than 60%, which highlights the importance of a correctly simulated surface albedo. The use of this coupled atmosphere–snow RCM offers new perspectives in the study of the Greenland surface mass balance due to the represented feedback between the surface climate and the surface albedo, which is the most sensitive parameter in energy-balance-based ablation calculations.  相似文献   

In Situ Measurement of Electro-Osmotic Fluxes and Conductivity Using Single Wellbore Tracer Tests     

Walt W. McNab Jr.  Roberto Ruiz 《Ground Water Monitoring & Remediation》2001,21(4):133-139

Electro-osmosis (EO), the movement of water through porous media in response to an electric field, offers a means for extracting contaminated ground water from fine-grained sediments, such as clays, that are not easily amenable to conventional pump-and-treat approaches. The EO-induced water flux is proportional to the voltage gradient in a manner analogous to the flux dependence on the hydraulic gradient under Darcy's law. The proportionality constant, the soil electro-osmotic conductivity or k_eo, is most easily measured in soil cores using bench-top tests, where flow is one-dimensional and interfering effects attributable to Darcy's law can be directly accounted for. In contrast, quantification of EO fluxes and k_eo in the field under deployment conditions can be difficult because electrodes are placed in ground water wells that may be screened across a heterogeneous mixture of lithologies. As a result, EO-induced water fluxes constitute an approximate radial flow system that is superimposed upon a Darcy flow regime through permeable pathways that may or may not be coupled with hydraulic head differences created by the EO-induced water fluxes. A single well comparative tracer test, which indirectly measures EO fluxes by comparing wellbore tracer dilution rates between background and EO-induced water fluxes, may provide a means for routinely quantifying the efficacy of EO systems in such settings. EO fluxes measured in field tests through this technique at a ground water contamination site were used to estimate a mean k_eo value through a semianalytic line source model of the electric field. The resulting estimate agrees well with values reported in the literature and with values obtained with bench-top tests conducted on a soil core collected in the test area.  相似文献